Logarithmic Correlations in Quenched Random Magnets and Polymers 
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It is argued that logarithmic factors multiplying power law behavior are to be expected in the 
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Logarithmic corrections to power law behavior at or 
near a critical point are normally associated with the 
presence of operators which are marginally irrelevant un- 
der the rcnormalisation group. The amplitudes of such 
logarithmic terms vanish if the critical theory happens 
to lie exactly at the fixed point. However, it has been 
pointed out [Q that the structure of two-dimensional con- 
formal field theory, which describes the fixed point in a 
wide range of isotropic systems with short-range interac- 
tions, allows so-called logarithmic operators in the fixed 
point theory itself, whose operator product expansions 
and correlation functions contain logarithms as well as 
pure power behavior. 

Such operators should not occur in unitary conformal 
field theories, such as correspond to pure critical sys- 
tems with positive Boltzmann weights, for the following 
reasons. Any real local density <i>(r) may always be ex- 
panded in a series of scaHng operators <i>(r) = 0^0; (r), 
so that its two-point correlation function in R*^ takes the 
form 



A. 



(1) 



where the scaling dimensions the eigenvalues of 

the dilatation operator T). The off-diagonal terms with 
Xi 7^ Xj vanish by conformal symmetry. In addition, the 
diagonal coefficients An are all non-negative as a conse- 
quence of reflection positivity. Thus, in such a theory, 
only pure powers, or sums of them, should ever appear 
[^. However, in non-unitary conformal field theories, it 
is possible for some of the An to be negative, in which 
case cancellations may occur in (Q). In particular, if two 
scaling dimensions Xi and xj become degenerate in such 
a way that An ^ —Ajj oo with Anlxi —Xj) remaining 
finite, the leading terms will cancel leaving a logarithmic 
term proportional to r~'^^' Inr. 

The main point of this paper is that such logarithmic 
terms should always occur in models which are described 
by non-trivial fixed point theories with a trivial parti- 
tion function Z — 1, in an arbitrary number of dimen- 
sions. These include all sorts of quenched random sys- 
tems including percolation, as well as problems such as 



self-avoiding random walks describing polymer statistics. 
Such logarithmic operators have already been identified 
in some systems jsj which may be analysed using super- 
symmetry to perform the quenched average, and there 
has been considerable work 1 1| on the general structure of 
such logarithmic CFTs. We shall use instead the replica 
formalism, which has the advantage of being more gener- 
ally applicable, as well as providing a simple way of seeing 
how the logarithmic terms appear in the limit. It will be 
seen that logarithms appear as the leading contributions 
in some quite simple two-point correlation functions, as 
well as completely generally as subleading terms in higher 
order ones. 

The analysis depends on some very simple properties 
of conformal field theory which are valid in all dimen- 
sions, not only d = 2. That is, they should hold for all 
critical theories of the above types as long as d is below 
the upper critical dimension. A crucial role is also played 
by the symmetry group of the theory, be it the permu- 
tation group Sn of the replicas, or 0{n) for polymers. 
Rather than work in complete generality we shall give 
two illustrative examples which are moreover physically 
relevant. 

Consider first a quenched random fcrromagnct, for 
example an Ising model with random bond disorder. 
Such systems may be described generically in the con- 
tinuum limit by a reduced hamiltonian H = H'^ + 
J m{r)E{r)d'^, where iJ" describes the pure system and 
m{r) is a quenched random variable coupled to the lo- 
cal energy density E{r). The disorder is supposed to 
have short-range correlations only, so that m{r) =0 and 
m{r)m{r') = gSrr'- Introducing replicas a — 1, . . . ,n in 
the standard way and averaging over the disorder gives 
the replica hamiltonian 



where the higher cumulants have been dropped on the ba- 
sis that they are irrelevant in the RG sense. The terms 
with a = b are discarded because they can be further re- 
duced using the operator product expansion (OPE) in the 
pure theory which takes the form Ea-Ea — l + boEa-\ . 
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The dimension of i? is = d — 2x'^, and it is relevant 
at the pure fixed point if d > 2a;^. This is equivalent 
to the well-known Harris criterion dv'^ < 2 If is 
small, it is possible to develop a perturbative RG in pow- 
ers of this variable and to locate a random fixed point to 
which critical systems with sufficiently weak randomness 
should flow ^,01 . The shift in the scaling dimension x^j, 
of some scaling operator </> from the pure to the random 
fixed point is given by M 



where b and 6^ are the coefficients in the OPEs $ • $ = 
1 + 6$ + ■ • • and (f)-^ — b^(j) + - ■ ■, and $ is the perturbing 
operator, in this case J2a^b ^aEt, which, using the above 
OPE of the pure theory, gives b = 4(n - 2) -I- 26§. 

Now consider the case where cj) is itself the replicated 
energy density Ea- The main point is that all these n 
operators are degenerate at the pure fixed point, and, to 
carry out the perturbative expansion, they should be de- 
composed into irreducible representations of the replica 
permutation group Sn- In this case this is straightfor- 
ward: the combination i? = E'a is a singlet, leaving 

Ea = Ea - (1/n) J2b ^b, with J2a ^o- = 0, transform- 
ing according to an (n — l)-dimensional representation. 
The OPE coefficients are easily found: bs = 2(n — 1) 
and b~ — —2, so that, close to n = 0, XE{n) ~ 
x% + \{l-n)yg + 0{yl) and x~{n) ^ x% + \yg + 0{yl). 
Note that the two scaling dimensions agree at n = 0. 
This is in fact true to all orders, and is a consequence 
of the fact that the partition function of any CFT 
on Si X Sd-i takes the form Z = where g is a 

modular parameter. Since Z = 1 at n = 0, all scaling 
operators must arrange themselves to have multiplicities 
proportional to n as n — > 0. Thus the singlet operator E 
becomes degenerate with the (n — 1) operators Ea- As- 
sume that the difference between the scaling dimensions 
vanishes linearly in n, as it does to first order. 

The two-point functions therefore take the form, for 

{E{r)E{Q)) = n {{E^E^) + (n - 1){E^E2)) 

- nA(n)r-2^^(") (2) 
{Ea{r)Eam = (1 - l/n) {{E^E^) - {E1E2)) 

~ (1 - l/n)I(n)r^^'^E("^ (3) 

where replica symmetry has been used in the first equali- 
ties. The amplitudes A{n) and A{n) have finite limits at 
n = in the pure theory and therefore should remain so 
at the random fixed point. For consistency between (|^) 
and (^) A{Q) — A{Q). Subtracting the two expressions 
and letting n — > 0, 



{E{r)){E{Q)) = \iu,{Ei{r)E2{r)) 



The left hand side of this equation gives the physi- 
cal meaning of the replica expression {E1E2): it is the 
quenched average of the disconnected part of the energy- 
energy correlation function. Note that the connected 
piece {EE) - (E) (E) = (EiEi) - {E1E2) is given by the 
n ^ limit of both equations above and contains no log- 
arithm. Although the amplitude A(0) is non-universal, 
the ratio of the quenched average of the disconnected 
piece to that of the connected piece is asymptotically 
proportional to Inr with a universal amplitude. 

Such a phenomenon does not occur for the correlation 
functions of the local magnetisation a ^ . Taking the ex- 
ample of an Ising model, the replicated hamiltonian is in- 
variant under independent Z2 transformations da —^ —Oa 
in each replica, as well as replica permutations. Under 
this larger symmetry, the operators {era} transform irre- 
ducibly, so that X(j(n) — xa{n) for all n. However, as ar- 
gued below, logarithms do appear in higher-point corre- 
lations of (J. Logarithms also appear in certain two-point 
functions of composite operators Ea^Ea2 ■ ■ -Ea^, which 
correspond to higher moments {EE)'p of the two-point 
function (and which exhibit multiscaling) . These arise 
]lOt from a similar mechanism to the above: the singlet 
operator {J^b^b)^ becomes degenerate with the {n — 1) 
independent operators {J2b^by^^i^a — / J2b ^b) , 
but only exactly at n = 0. Note also that it has been 
observed |jri| that logarithmic operators whose 2-point 
functions behave like {D{r)D{0)) r~^^ (In r + const.) 
are always accompanied by another operator C such that 
(C(r)L>(0)) - and (C(r)C(O)) = 0. In this exam- 
ple, C = E and D ~ Ea, in the limit n = 0, so that D is 
a reducible operator under the replica symmetry. 

As a second example, consider self-avoiding walks, 
which are a model for polymer statistics, and are de- 
scribed by the n ~> limit of an n-component 
field theory with a hamiltonian with 0{n) symmetry 
/ (Ea((V0a)2 + m202) + g ^2^, <j>l4>l) d^r . Consider 
the composite operators of the form : 0a 0b ■ (the nor- 
mal ordering subtracts the disconnected part {4>a4>b))- 
These may be decomposed under 0{n) into a singlet 
E = : 0^ : and a symmetric traceless part Eab ='■ 
(f>a(t)b ■■ -{^/n)SabT,c (t'c with {n - l){n + 2)/2 inde- 
pendent components. In general, these two irreducible 
sets carry different scaling dimensions, as can checked in 
perturbation theory in g: at one loop, the relevant OPE 
coefficients are bs = 8 + 4n and b~ = 8. Once again, they 



become degenerate at n = so that Z = 1 |13|. Their 
two-point functions have the form 



{E{r)E{0)) = 2nA(n)r'2a;E(«) 

{EabEcd{0)) = {SacSbd + hcSad ^ fjabScd.)A{n)r 



-2x~(n) 



2A{Q)x'e{Q) lnr/r-2^^(°) 



where the index structure is fixed by group the- 
ory and the amplitudes A{n) and A(n) are finite 
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at n = 0. Now consider the correlation function 
(: 4>a{r)(l)b{r) :: (/)a(0)(/)f,(0) :) and express this in terms 
of {EE) and {EE). There is then a singular term pro- 
portional to ((5ab(5cd/n)(^(n)r^2a;E(n) _ A{n)r~'^'''B^''^), 
which gives rise to a logarithm at n = 0. A similar 
argument may be made away from the critical point, 
where the two-point function of each scaling operator 
should scale as ^^^^'^"■'/^('"/C; where ^ is the correla- 
tion length. This implies that the generalized susceptibil- 
ities XE = Y.r{E{ri)E(r2)) and = ('^2)) 
should behave as Ci{n)£_'^^'^^'^"-\ where i £ {E,E} and 
Ce{0) = C~(0). 

The physical interpretation of this is as follows. The 
term proportional to Sab^ai in the full correlation func- 
tion counts all pairs of mutually avoiding self-avoiding 
loops (polygons on a lattice) one of which goes through 
ri, and the other through r2. The normal ordering sub- 
tracts off the disconnected piece which corresponds to the 
product of the independent sums over such pairs of loops, 
in which they are allowed to intersect. Therefore the re- 
maining connected part, to which the above is the leading 
contribution, corresponds to all configurations where the 
two loops intersect at least once. Denoting by cat the 
number of such pairs of loops, of total length N , per lat- 
tice site, the mapping to the 0{n) model is such that 
its generating function N'^cnV^ is proportional to 
lim„_o(xB - Xg)/", with ^ ex (j/c - v)'" ■ This implies 
that CAT - N^d-'^''''>-^n^hiN = iV"-V^lniV, where 
ji = y^^ . For monodisperse loops, there is an additional 
factor of N~^. This agrees with an alternative method: 
the free energy per site f{n) of the 0{n) model scales as 
B{n){xc — a;)^""*^"-', where i?(0) = 0. Since each loop ap- 
pears with fugacity n, the second derivative df /dn\n=Q 
counts pairs of intersecting loops, and its leading singu- 
larity |l| is of the form B' {0)a' {0){xc - \ii{xc - 
x), consistent with the above result. 

Although we have given two examples of how loga- 
rithmic factors occur in two-point functions, it is worth 
pointing out that they are ubiquitous in four- and higher- 
point functions in fixed point theories with Z = 1. This 
is because conformal invariance fixes the coefficient of the 
stress tensor T in the OPE of an operator (j) with itself 
to have the form (j) ■ (j) = 0^(1 -I- {x^/c)T + ■ ■ •), where all 
the co-ordinate dependence and the complicated index 
structure has been suppressed [p^ . Here defines the 
normalisation of and c is the central charge, defined in 
arbitrary dimension in terms of the two-point function 
(r(r)r(O)) cx c/r"^, where again the index structure has 
been suppressed. Substituting these terms of the OPE 
into the 4-point fimction {4>{'''i)4'{i'2)4'{i'3)4'i'''i)) in the 
limit where rj = ri2r^4^/ri^r2A ^ 1 gives a term pro- 
portional to (ri2r34)~^^* (a^/c)ry''. In a theory such as 
those discussed above, with c = 0, this raises a poten- 
tial difficulty The resolution is that, in general, 
there are other operators whose scaling dimensions be- 



come degenerate with that of T as n 0. For the 
quenched random ferromagnet, these are the operators 
Ta ^ Ta ~ {^/n) X^b^ftj constructed in a similar way to 
the Ea above. For the 0{n) model, they are traceless 
part of Tj^^ = 9^(/)a9'^0b + • ■ ■, whose trace over a = b 
gives the stress tensor T'^". These operators give rise to 
logarithmic terms of the form 77"^ In rj in the four-point 
functions at n = 0. 

However, the paradox may be resolved in another way. 
The physical correlation functions may be such that 
— > as c — > 0. This in fact happens in the g-state 
Potts model where all connected correlations vanish lin- 
early at g = 1, so that oc c. The connectivities of the 
percolation problem are given by the derivatives with re- 
spect to (? at q = 1, and are finite. But in this case there 
are no logarithmic terms of the above form 

It would be very interesting to understandjiow the 
(n — 1) bosonic partners of the stress tensor like T become 
its fermionic partners in those cases where the quenched 
average can be performed using supersymmetry [T^ . 

We have argued that logarithmic factors multiplying 
power law singularities are ubiquitous in critical systems 
described by effective theories with Z ~ 1. In general 
they hide themselves in non-leading corrections to higher- 
point correlations, but there are examples where they 
appear in two-point correlations and susceptibilities. Al- 
though all conclusions are for general d, it is possible to 
check some of them in d = 2 for percolation and the 
0{n) model [ [l9| ], where many of the critical four-point 
functions are known exactly. 
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